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ABSTRACT 

Mean-motion resonances between a Keplerian disc and an orbiting companion are anal- 
ysed within a Hamiltonian formulation using complex canonical Poincare variables, 
which are ideally suited to the description of eccentricity and inclination dynamics. 
Irreversibility is introduced by allowing for dissipation within the disc. A method is 
given for determining the rates of change of eccentricity and inclination variables of 
the disc and companion associated with resonances of various orders, including both 
reversible and irreversible effects, which extend and generalize previous results. Pre- 
liminary applications to protoplanetary systems and close binary stars are discussed. 
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1 INTRODUCTION 

1.1 Satellite— disc interactions 

The gravitational interaction of a gaseous or particulate disc with an orbiting companion is of general interest in astrophysics. 
Important examples include protoplanetary systems involving one or more planets orbiting within a gaseous disc, accretion 
discs in close binary stars, and systems of planetary rings and moons. 

As in the planetary theory of classical celestial mechanics, the disc and companion experience deviations from perfect 
Keplerian motion around the central mass as a result of their mutual gravitational perturbations. Sufficiently close to the 
companion's orbit, the behaviour of the disc can be complicated and highly nonlinear, leading in many cases to an exclusion 
of the disc from the coorbital region. The extent of such a gap depends on the ratio of the mass of the companion to that of 
the central object, and on the properties of the disc. Away from this region, the interaction can be analysed by perturbation 
methods and has two principal forms. 

Secular interactions (of first order in the perturbation) arise from the mass distribution of the companion averaged around 
its orbit. These contribute to the precessional behaviour of the system and allow slow, reversible exchanges of eccentricity 
and inclination between the companion and the disc. Such behaviour is analogous to the Laplace-Lagrange secular theory 
of planetary systems in which, for example, Jupiter and Saturn undergo large -amplitude oscillatory exchanges of eccentricity 
and inclination over tens of thousands of years (e.g. Murray fc Dermott 19991 ). 

Mean-motion resonances depend specifically on the periodic nature of the companion's orbit. The fluctuating forces 
induced by the companion generally give rise to rapid oscillations of small amplitude in the disc that are of little consequence. 
However, where a commensurability occurs between the orbital periods of the companion and a particle in the disc, a strong 
localized interaction occurs. Through the intervention of collective effects and dissipation in the disc, a different kind of secular 
behaviour is induced at second order in the perturbation. This leads to irreversible changes in semimajor axis, eccentricity 
and inclination of the disc and companion, rather than just precessional behaviour. 



1.2 Corotation and Lindblad resonances 



Mean-motion resonances in sat ellite-disc interactions were studied in de tail in an influential series of papers by Goldreich & 
Tremaine. Early in this series (|Goldreich fc Tremainel Il978l . Il979l . Il980l ) the authors considered a two-dimensional problem 
in which a circular disc is perturbed by a companion with a slightly eccentric orbit. To calculate the linear response of the 
disc, they decomposed the perturbing potential by Fourier analysis into a series of rigidly rotating components with various 
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amplitudes, azimuthal wavenumbers and angular pattern speeds. In a system with Keplerian orbits, the angular pattern 
speeds that occur are rational multiples p/q of the mean motion of the companion, the amplitude scaling with the \p — q\th 
power of the eccentricity. 

The linearized equations governing the response of the disc at a given azimuthal wavenumber and angular pattern speed 
allow for the propagation of free waves in certain intervals of radius. Three important radii are the corotation radius, where 
the wave frequency as measured in a frame rotating locally with the disc vanishes, and the two Lindblad radii, where the same 
quantity equals (plus or minus) the epicyclic frequency. In the simplest case of a two-dimensional non-self-gravitating disc, 
only one wave mode is permitted. This 'density wave' combines inertial (epicyclic) and acoustic behaviour and can propagate 
only where the wave frequency seen by the disc exceeds the epicyclic frequency in magnitude. It therefore propagates interior 
to the inner Lindblad radius and exterior to the outer Lindblad radius, but is evanescent between the Lindblad radii where 
the corotation radius lies. 

The forcing of disturbances by a perturbing potential component is effectively localized to the neighbourhoods of the 
corotation and Lindblad radii. Beyond the Lindblad radii, the forcing is ineffective because of the very limited overlap between 
the potential and the free waves, which oscillate rapidly in radius. Despite the evanescent nature of the free waves near the 
corotation radius, it plays an important role as a singularity of the linearized equations. The forcing can then be said to occur 
at corotation and Lindblad resonances, which can be identified with the mean-motion resonances of celestial mechanics. The 
reason for this is that, in a Keplerian system where the epicyclic frequency equals the orbital frequency, any radius where 
the orbital frequency is commensurate with the mean motion of the companion serves, in principle, as both corotation and 
Lindblad resonances for some potential components. 

At a corotation resonance the linearized response of the disc at first order in the perturbation is singular but can be 
resolved by including a viscosity. At second order in the perturbation, the companion exerts a secular torque on the disc as it 
generates a localized, non-wavelike disturbance that transfers angular momentum steadily to the disc through viscous stresses. 

At a Lindblad resonance the localized response is regularized by collective effects such as pressure, self-gravity or viscosity, 
and is generally in the form of an attenuated wave propagating away from the resonance. At second order in the perturbation, 
and provided that the w ave is dissipated somehow with in the disc, the companion again exerts a secular torque in launching 



the disturbance (see also iMever-Vernet fc Sicardv 



In a three-dimensional disc, not considered by 



19871) 



Goldreich fc Tremaine! ( 19791 ). a vertical resonance is similar to a Lindblad 



resonance but involves forcing and motion perpendicular to the plane of the disc, so that a bending wave is excited rather 
than a density wave. 

Convenient formulae are available for the corotation and Lindblad torques, which have an important role in determining 
the orbital evolution of the companion and the surface density distribution of the disc. The subseq uent literature has made 
extensive use of these expres sions, as well as providing some modifications to the torque formulae (|Wardlll989l : lArtvmowiczl 
19931 : iQgilvie fc LuboJ2003l ) . Together, these results provide the basis of our understanding of planet-disc and other satellite- 



disc interactions. 



1.3 Alternative treatments of mean- motion resonances 

It is perhaps not widely appreciated that the approach described above is inadequate for general satellite-disc interactions. A 
knowledge of the resonant torque (by which is meant the component of the torque perpendicular to the disc) determines one 
component of the rates of change of the orbital angular momentum vectors of the disc and companion and also, by virtue of 
Jacobi's theorem, the rate of change of orbital energy of the companion; however, this information does not generally suffice 
to determine the secular rates of change of the relevant orbital elements of the disc and companion. Furthermore, the method 
does not apply to the general case in which the disc is eccentric, because the forced wave equations in an eccentric disc have 

not been analysed. 

Later in their series of papers, Goldreich fc Tremaine! 1 198lh and Borderies. Goldreich fc Tremaine] 1 19841 ) employed meth- 



ods of celestial mechanics to treat cases in which either the companion or the disc has an eccentric and/or inclined orbit. Using 
'simple artifices' they derived formulae for the rates of change of the canonical action variables associated with mean-motion 
resonances, in terms of the classical disturbing function, while avoiding an explicit discussion of collective effects. In principle, 
these expressions are considerably more powerful and general than the earlier torque formulae, yet they appear not to have 
been used in the subsequent literature. We have tried to make use of them in connection with eccentric planet-disc interactions 
and similar problems but found that they generally give results that are not of the desired form. 

One aspect of their method is that it isolates, in the traditional way, terms of different angular arguments in the disturbing 
function. An implicit averaging is therefore involved over the mutual apsidal and nodal precession of the disc and companion. 
At least in some contexts, this procedure is inappropriate because the mutual precession may occur on a timescale comparable 
to that of the secular evolution induced by mean-motion resonance s, or indeed may not occur at all if the disc and companion 



are engaged in a secular normal mode fe.g. iLubow fc Ogil vie 2001). The formulae therefore describe the rates of change of a 
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e and i averaged over the timescales of mutual precession, while what are generally desired are the rates of change of a, e, i, 
zu and 0. averaged only over the relative mean motion^ 

A second aspect is the way that dissipation is introduced into the analysis, which is essential to capture the irreversible 
natur e of the effects of mean-motion resonances. This is done either by including ad hoc damping terms (jGoldreich fc Tremaine 



1981) o r is avoided altogether by using Landau 's prescription in which the disturbing function is given a slow exponential 



growth ijBorderies. Goldreich fc Tremaindll984h . Although mathematically convenient, these approaches can give misleading 
results. The effects of mean-motion resonances depend on the fact that the disc conserves angular momentum while it dissipates 
energy, and the damping terms added to the equations ought to respect these properties. The Landau approach is problematic 
because the first-order oscillatory disturbances in the disc are permitted to grow exponentially rather than being steadily 
dissipated, making it difficult to calculate the second-order secular changes correctly. 

The purpose of the present paper is to revisit the analysis of mean-motion resonances in satellite-disc interactions and 
to derive useful and general formulae for the associated evolution of orbital elements of the disc and companion. We employ 
methods of celestial mechanics and treat the collective effects of the disc in an approximate way that is adequate for our 
present purposes. The remainder of the paper is organized as follows. In Section [2] we consider the restricted problem in which 
the companion has a prescribed orbit. We develop the theory using a Hamiltonian formulation of the equations, and describe 
the role of collective effects. In Section[3]we consider the coupled satellite-disc system consistently and determine the feedback 
of mean- motion resonances on the companion. Preliminary applications are worked out in Section [4] and a summary of the 
main ideas and results follows in Section [5] 



2 THE RESTRICTED PROBLEM 
2.1 Hamiltonian formulation 

We consider a thin disc around a central mass Al, and initially neglect any collective effects. The disc is perturbed by an 
orbiting companion of mass M' . In the restricted problem the perturber has a prescribed periodic orbit that is unaffected by 
the disc. This approximation is usually applicable, for example, to accretion discs in binary stars. In this case the influence of 
the companion on a test particle in the disc can be described through a disturbing function that is a specified function of the 
position of the particle and is also a periodic function of time. 

In order t o treat the dynam ical equations in the most efficient way, we adopt a Hamiltonian formulation using canonical 
variables (e.g. Morbidelil 2002 ). The mean longitude A is conjugate to the energy-related action variable A = (GMa) 1 ^ 2 , 



where a is the semimajor axis. The other two generalized coordinates and conjugate momenta, which relate to eccentricity 
and inclination and need not be in action-angle form, are written as q a and p a , with a — {1,2}. Hamilton's equations for a 
test particle are 

; dH : dH . dH . dH 

A =5A' A = "9A' qa = dp~ a > Pq = ^' (1) 
where H is the Hamiltonian per unit mass[f] The unperturbed (Keplerian) Hamiltonian 
(GM) 2 GM 

2A2 - = 2a (2) 
gives rise to a mean motion 

• ... dH {GM) 2 /GM\^ 2 

A = n(A) = — = L—L = I— j . (3) 

The disturbing function is periodic in the mean longitudes A and A' of the test particle and the perturber, and can 
therefore be expanded in a double Fourier series. We consider a Hamiltonian of the form 

H = Ho ( A) - e Re [R( A, q a , p a ) e i0 ] , (4) 

where e = 1 is a bookkeeping parameter used to identify effects of various orders in the perturbation, R is a complex potential 
amplitude (proportional to M') and <f> — kX + k'X' is a potentially resonant angle with integer coefficients k and k □ In the 
restricted problem A' = n't + e' is a prescribed linear function of time. The double Fourier series in fact gives rise to an 
infinite number of perturbing terms of this form, which might be labelled using the notation Rk t y an d 4>h,h' ■ We may consider 



1 We use the standard notation for the Keplerian orbital elements (e.g. iMurrav fc D crmott 1999]): a, e, i, vj and O denote the semimajor 
axis, eccentricity, inclination, longitude of pericentre and longitude of ascending node, while A is the mean longitude. The secular 
contribution to the evolution of A is of less interest here. 

2 The true Hamiltonian of the system is J H dm, where dm is a mass element of the disc. 

3 Note that it is permissible to replace R with R* if the signs of k and k 1 are both changed. 
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them separately because their second-order secular effects will sum in quadrature owing to the phase relations between them. 
Hamilton's equations are then 



A = n(A) -eRe 



dR i, 

P 

OA 



A = eRe [ifeiJe 10 ] , 



q a — —eRe 



dR j 4 



p a = eRe 



dR 
dq a 



(5) 
(6) 

(7) 
(8) 



An important point is that the apsidal and nodal precessional longitudes vj and O appear in R (through its dependence 
on q a and p a ) and not in the exponential. Terms in the disturbing function with the same values of k and k! but different 
dependences on vo and Q are treated together, not separately. To isolate such terms would be appropriate only if we were to 
average the behaviour of the system over the mutual apsidal and nodal precession of the disc and companion. 



2.2 Perturbative expansion to second order 

We expand the dynamical variables in the form A = Ao + eAi + e 2 A2 + etc., with the intention that Ao describes 
the unperturbed disc, Ai the periodic oscillation induced by the perturbation at first order, and A2 contains second-order 
variations. The effects that we are looking for will appear as secular changes of A2, etc0 
At leading order all variables are independent of time except for 



A = n(A ) = n , 



(9) 



which describes the unperturbed mean motion of the disc. It is convenient to use either the unperturbed semimajor axis ao 
or the related action variable Ao = (GMao) 1 ^ 2 to label the orbits in the disc. Then 



bo = kno + k n 



(10) 



is a constant depending on Ao. By definition, it vanishes at the location of the mean-motion resonance where no/n = —k'/k 
is a certain positive rational number. 

At first order we find from Hamilton's equations 



dA 



Re 

Ai = Re [ikRe^ ] 



dR j 
8 A 



00 



q a i = —Re 



OR 

dp a 



,i0O 



p a i = Re 



dR 
dq a 



,i0o 



(11) 
(12) 

(13) 
(14) 



where all terms on the right-hand sides are evaluated on the unperturbed solution, i.e. at A = Ao, q a = qao, Pa = PaO- The 
solution of equations (|12[l - (|14p is 



Ai = Re [fkRe^ ] 



q a i 



Pal 



Re 



if 



Re 



dR 

dpa " 

dR 



.100 



-if 



dq a 



100 



(15) 
(16) 

(17) 



where 



(18) 



4 This form of perturbative expansion is valid only for a limited time interval, because A2 will eventually become comparable to Ao. 
A more sophisticated approach would be to use the method of multiple timescales, which would yield a uniformly asymptotic solution. 
However, the notationally simpler approach taken here allows us to calculate the secular rates of change correctly. 
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Equation then becomes 
Ai = Re 



dA J OA 



Re 



for which the solution is 



Ai = Re 



(19) 



(20) 



and then we find (j>i = k\i. To this order, therefore, the variables undergo forced oscillations around their leading-order values 
with angular frequency (j>o- Four arbitrary constants could be included in the solution, corresponding to small shifts in the 
mean values of the variables, but these are best absorbed into the leading-order terms. The amplitude of the oscillations 
diverges at the resonance where (fio = 0. 

At second order the equation for A2 is not required, but we have also 



A 2 =Re 



ifc 



dR , OR OR ., n 

dA Al + d^ q01+ d^ P!31+1 * lR 



q a 2 



-Re 



Pa2 = Re 



2 R 
dp a d A 

d 2 R 



Ai + 



dq a dA 



Ax + 



d 2 R 

dp a dqp 

d 2 R 



d 2 R . ± dR 

dpadpp dp a 



i©0 



dq a dq, 



-qgi 



d 2 R 

dq a dpp 



Ppl + 101 



dR 

dq a 



,i0o 



(21) 



(22) 



(23) 



where, again, terms on the right-hand sides are evaluated on the unperturbed solution, and we adopt the convention that 
a summation over (3 — {1,2} will be implied wherever it appears. The right-hand sides, involving products of oscillating 
quantities, contain both non-oscillatory terms of zero frequency and oscillatory terms of frequency 2</>o- We are interested in 
the secular rates of change of A, q a and p a at second order in e, which we denote by angle brackets. Making use of the relation 



Re [A e i(p0 Re [B e i0 °] ] = ^Re [AB e 2 ^ + A* B] 



and discarding the oscillatory terms, we find 



(A) = iRe 



dR* dR dR* dR 



{q a ) = -- Re 



if 



dqp dp/3 
d 2 R* 



dpp dqp 
dR d 2 R 



i^(/l*! 2 ) 



dR 



Re 



if 



dp a dqp dp/3 

d 2 R* dR 



dp a dp dqp 
d 2 R* dR 



+ k dA 



fR 



dR* 



dq a dqp dp dqadpfj dqp 



+ k 



d_ 

dA 



fR 



dp, 

dR* 
dq a 



(24) 



(25) 



(26) 



(27) 



So far we have not used the fact that / is a real quantity; the reason for this will become apparent later. Away from the 
resonance, / is real and finite and equations (|25l) - (|27|l become 

dU 



(A> = 
(q a ) 



&H 

dp a 



= 0, 



(Pa) = -75—. 

aq a 



where 



-if 



{ dR* dR 
I dqp dpp 



dR* dR 
dpp dqp 



k^if\R?) 



(28) 
(29) 
(30) 

(31) 



is a real Hamiltonian. 

The implication of our analysis is that each oscillatory term Rk,k' m t ne disturbing function gives rise to an oscillatory 
behaviour in the disc at first order and to a secular behaviour at second order. Since it is described by a Hamiltonian H.h,y 
that is independent of A, this secular behaviour concerns only the eccentricity and inclination variables and is reversible and 
precessional in character. It is of interest mainly because it diverges in the vicinity of mean-motion resonances where / — > 00. 
The resolution of this singularity, which results in non-Hamiltonian, irreversible secular behaviour, is discussed in Section \2. 51 

The above procedure does not work for the secular part f?o,o of the disturbing function, for which cj>o vanishes identically. 
However, this term contributes directly to the secular dynamics at first order in the perturbation, giving rise to the familiar 
precessional effects, and we are not concerned with its smaller, second-order counterpart. Second-order effects are of interest 
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only because of their amplification by mean-motion resonances. The secular dynamics is therefore adequately described by a 
Hamiltonian consisting of —Ro,o together with the sum of Ttk,k' f° r an Y important resonances {k, k'}, subject to the resolution 
of the corresponding singularities. 



2.3 Complex canonical variables 

Although Hamilton's equations employ real variables, complex variables are natural in c elestial mechan ics, especially for 
describing precessional behaviour. The idea of complex canonical variables was employed by Strocchi ( 19661 ) in order to unify 
the descriptions of classical and quantum-mechanical systems. Suppose we have a Hamiltonian system with canonical variables 
q a and p a and Hamiltonian H (? a ,p Q )EI Then (q a ,p a ) can be replaced by (z a , z„), where 



z a = — ( 9o 



+ VPa) 



(32) 



(In order for the dimensions to match, q a and p a should be variables of 'rectangular' rather than 'action-angle' form 
Hamilton's equations then combine in the compact form 

dH 



dz* 



(33) 



in which H is now regarded as a real- valued analytic function of z a and z a , treated as algebraically independent variables. 
(The conjugate equation — ii* = dH/dz a need not be considered separately, as it follows directly from the complex conjugate 
of equation (|33[) ; z a and z* a are 'conjugate' variables in both senses of the word.) 

The equations derived so far in this paper are valid for any choice of canonical eccentricity and inclination variables q a 
and p a used in conjunction with A and A, such as modified Delaunay variables (e.g. Morbidelli 20021 ) . Indeed the quantity 

dR* dR OR* dR 



= {R*,R} 



(34) 



dq dpp dpp dqp 

that appears in equation (|31l) is a Poisson bracket and is therefore invariant under canonical transformations. (The variables 
A and A do not appear in this Poisson bracket because ii is independent of A.) We will use the canonical variables 

(35) 
(36) 
(37) 
(38) 

which are identical to the rectangular variables of Poincare] (1_892), but rotated through an angle of ty/2. When e <C 1 and 



9i = 


{2A[1- 


-(1-e 2 ; 


1/21 \V2 
' | COS?a7, 


Pi = 


{2A[1- 


-(1-e 2 


1/21 1 I/ 2 . 

' | sinter, 


<72 = 


[2A(1 - 


e 2 ) 1/2 (l 


— cosi)] ^ cosfi, 


P2 = 


[2A(1 - 


e 2 ) 1/2 (l 


\1 1 / 2 

— cosi) sinS7, 



i 1, these variables reduce to (qi,pi) ~ A 1 / 2 (ecosro, esinro) and (172, P2) ~ A 1//2 (i cos O, i sin fl). showing their relation to 
the Cartesian components of the eccentricity and inclination vectors (e.g. iMurrav fc Dermottlll999l h 
Strocchi's transformation leads us to consider 
, 1/2 



zi = {A [l 

Z2 = [A(l - 



2\1/21 1 1 

e ) J) 



(39) 



(1 

e 2 ) 1/2 (l-cosi)] i/ % in , (40) 
which may be called complex canonical Poincare variables and were also employed bv lLaskar fc Robutel (1995)0 Note that 



L = A - \ Zl \ 2 = [GMa(l - e 2 )] 



1/2 



is the total specific angular momentum and 



L z = A — |zi| 2 — \z 2 \ z = [GMa(l - e 2 )] 



1/2, 



(41) 



(42) 



is its component perpendicular to the reference plane. Since A is the specific angular momentum of a circular orbit of radius 
a in the reference plane, \zi\ 2 + \z2\ 2 is the specific angular momentum deficit or AMD ( Laskarlll997l ). It consists of two parts, 



\zi\ 2 associated with eccentricity and \z 2 \ 2 with inclination, and is a positive-definite measure of the deviation of an orbit from 
circularity and coplanarity (with respect to the arbitrary reference plane). 

We now consider the complex potential amplitude R as an analytic function of A, z a and z a satisfying the analytic 
properties 



In this paragraph we assume that a labels all the degrees of freedom, not just a subset. 

They can alternatively be obtained from the modified Delaunay variables via the transformation z a 



1/2, 
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dR* t dR \ * _ dR* 

i ) dZa ' V dZa J dz a ' 

The Poisson bracket of equation (|34|l becomes 
. f dR* dR dR* dR ' 



{R*,R} 



\ dzp dz*^ dz*p dzp 
and so equation (|31[) for the second-order Hamiltonian reads 



f 



dR 



dzp 



dR 



dz* 



(43) 



(44) 



(45) 



It will be seen below that the terms in the first bracket, which arise from the dynamics of the eccentricity and inclination 
variables z a , correspond in some sense to Lindblad and vertical resonances, while the term involving a derivative with respect to 
A, which arises from the dynamics of the mean-moti on variables A and A, corresponds i n som e sense to corotation resonances. 
A similar distinction is apparent in the analysis of iBorderies. Goldreich fc Tremaind (|l984T ). We therefore refer to the first 
and second parts as the Lindblad /vertical term and the corotation term, respectively. 



2.4 The disturbing function in complex canonical Poincare variables 



The disturbing function is usually expanded in the variables e and s = sin(i/2). In converting these into complex Poincare 
variables it is useful to define the dimensionless complex eccentricity and inclination variables 

/ 2 \ 1 / 2 

£=[j) zi, (46) 



-(if- 

and then to note that 



2\l/2 



2se =1(1- ||£| 



-1/2 



(47) 

(48) 
(49) 



Therefore the magnitudes of £ and I are directly related to the eccentricity and inclination and are approximately equal to 
them when small, while the phases of £ and I are the corresponding precessional longitudes zu and fi. 

Expansions for R in terms of £ and T up to fourth d egree in eccentricity and inclination are given for the most important 
commensurabilities in Appendix |X] Lask ar fc Robutell (|1995l ) provide a method for generating such expansions, and prove 
theorems regardi ng their general alg e braic form. In fact, the expressions in Appendix [X] were obtained simply by rewriting 
the expansion of Murray fc Dermott (1999) in terms of £ and X. (Note, however, that primed variables always refer here to 
the companion, regardless of whether the resonance lies interior or exterior to it.) 

For the purposes of illustration, consider a first-order interior j : j — 1 resonance with j > 2. By this notation we mean 
that the resonant orbit in the disc is interior to the perturber and has a mean motion j / (j — 1) times that of the companion. 
We define the ratio of semimajor axes a — a/ a' — [(j 
the sum of direct and indirect terms in the form 

GM' 



l)/i] 2/3 . [Murray fc DermottJ |l999h give the disturbing function as 



K 



(TZd + aTZ E ) ■ 



(50) 



The terms of lowest degree in eccentricity and inclination for this commensurability are (jMurrav fc Dermottl fl999. Tables B.4 
and B.5) 



Tin = e/27 cos[jA' - (j - 1)A - w] + e'/3i cos[jA' - (j - 1)A - a/], 
TZe = — 2e'cos[jA' — (j — 1)A — izr']8j t 2, 
where 8 is the Kronecker delta and 



d 



(51) 
(52) 

(53) 
(54) 



are given in terms of Laplace coefficients. Comparing with equation Q, we may therefore identify the integer coefficients 
k = j — 1 and k' = —j, the corresponding frequency <j)o = (j — l)n — jn' and the complex potential amplitude 
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R; 



GM' 



[faS + (/si - 2q5 3 , 2 )£:'] 



(55) 



The equivalent expression correct to fourth degree in eccentricity and inclination is given in equation (|A3f) . which compactly 
combines all terms in Tables B.4 and B.5 of lMurrav fc Dermottf (|l999h . Note that the inclination terms vanish when 1=1', 
i.e. when the disc and perturber are coplanar. 



2.5 Resolution of the resonant singularity 



As a simple but important example, at a second-order interior resonance with a perturber having a circular and coplanar 
orbit, such as the 3 : 1 resonance that occurs in accretion discs in close binary stars of sufficiently small mass ratio, the 
dominant term in the disturbing function has R oc £ 2 for e < 1 and is therefore of the form 



R = g(\)zl 



0(z% 



(56) 



where, in fact, g is real (cf. equation lA5|) . According to equation (|45[) . the second-order Hamiltonian arising from this disturbing 
function is 



H = . 



X 2i i2 

9 - 



Q(z 4 



which gives rise to the dynamics 

dH 

dz* 



i<ii) 



bo 1 g 2 z 1 + 0(z 3 ). 



(57) 



(58) 



Since arg(zi) = m, this behaviour corresponds simply to apsidal precession at a rate —g 2 /cj>o, which is retrograde interior 
to the resonance (where <f>o > 0), prograde exterior t o the resonance (whe re 4>o < 0) and diverges at the resonance itself. 
This behaviour can be seen, for example, in Fig. 1 of iKnezevic et al.l (| 199 11 ). which shows how the mean apsidal and nodal 
precession rates of asteroids in nearly circular orbits, calculated to second order in the planetary masses, are affected in the 
vicinity of mean-motion resonances with Jupiter. The 3 : 1 and 5 : 3 resonances display the characteristic behaviour described 
here for second-order interior resonances. 

This divergence arises because the expansion we have adopted breaks down in the vicinity of the resonance. The behaviour 
of non-interacting particles and of continuous discs near a mean- motion resonance is different. Single particles avoid the 
divergence by undergoing libration, which means that the resonant angle, instead of circulating arbitrarily slowly, oscillates 
about an equilibrium value. Fluid or collisional particle discs, in which the intersection of streamli nes is resisted, can avoid 



the di vergence through the intervention of collective effects such as pressure, self-gravity or viscosity. Me yer- Vernet fc Sicardv 



1 1987) showed how any of these collective effects, in the case of a Lindblad resonance, acts to displace the resonant pole 
slightly away from the real axis, after a spatial Fourier transform of the linearized perturbation equations is carried out. We 
initially adopt a similar prescription, replacing equation (|18[) with 

1 



/ = 



(59) 



where s is a positive parameter with the dimensions of frequency. This replacement can also be motivated by Landau's 
approach, in wh ich the periodic disturbing function is 't urned on' slowly by including an additional time-dependence propor- 
tional to e s * (cf. Borderies. Goldreich fc Tremainei 1984 ). However, Landau's approach must be used with care and we return 
to address this issue more thoroughly in Section [2.61 

When / is treated throughout the analysis of Section [22] as a complex quantity, we obtain secular evolutionary equations 
at second order that are not of Hamiltonian form but allow for the irreversibility that is associated with the resonance. In 
particular, 



(A) 



( 


dR 


2 


dR 






dzfi 




dz* 





(60) 



1 



d 2 R* dR d 2 R* OR 



d 2 R dR* 



d 2 R OR* 



(fR dR * | f*R* dR 
4 dA \ dzSi dz*. 



4 J ydzZdzp dzfs dzadzp dz*p J 4 ydzadzp dz^ dzadz^ dzp 
where / = Im[/]; these expressions agree with dTL/d\ and dTL/dz*, if /i = 0. In the present example we obtain, for e< 1, 
6o + is 



(61) 



i(£0 = fg zi 



2 

g zi, 



(62) 



instead of equation (|58|l . Resolution of the resonant singularity therefore leads to a finite precession rate of —g 2 4>o/(<j>Q + s 2 ) 
together with a qualitatively different effect: an eccentricity growth rate of g 2 s/(<p 2 + s 2 ). The growth rate peaks at the 
resonance, the height and width of the peak depending on the parameter s, while the integrated growth rate is independent 



Satellite-disc interactions 9 



of s. In fact the growth rate can be represented as ivg 2 l 
broadened. We note also that 



*(0O) = 



dA 



5(A - A ) 



da 



S(ao — ao), 



as s — > 0, i.e. in the case of a resonance that is not significantly 



(63) 



where Ao (or &o) locates the resonant orbit. Unlike the Hamiltonian precessional behaviour, the eccentricity growth is irre- 
versible and corresponds to a growth of the AMD of the disc: 



dt 



2s 



2 2 

g \zi\ 



r, 2| |2 

2ng \zi\ 



(64) 



The physical resolution of the divergent precession rate is therefore a finite precession rate together with an irreversible 
growth or decay (in this case a growth of eccentricity). This behaviour is characteristic of an eccentric Lindblad resonance 
as described by iLubowl (|l99ll ). For the 3 : 1 resonance in a system with mass ratio q = M' /M < 1 we obta in g ~ 
1.727 q(GM/a' 3 ) 1/2 an d an eccentricity growth rate of 2.082 q 2 (GM/a' 3 ) 1/2 a 5(a - a ). This agrees with the result of lLubowl 



( 199lh as interpreted bv lGoodchild fc Ogilvie] (|2006h . '. iowever, neither of these papers considered the second-order precessional 

effect. 

A pplying equati on (14) of Borderies. Goldreich fc Tremainel 1 1984 ). however, gives exactly twice the eccentricity growth 
rate of lLubowl 1 1991 ). The reason for this is that their (modified Delaunay) variable F is proportional to e 2 for e < 1. In their 
analysis e 2 grows partly because the non-oscillatory part of e grows secularly, which is physically correct, but also because the 
oscillatory part of e grows proportionally to e st . In reality the oscillatory part of e has a steady amplitude owing to dissipation 
in the disc. Therefore application of Landau's prescription in this case gives a misleading result for the second-order behaviour 
of the disc. 



2.6 Collective effects in the disc 



We now attempt to give a justification for the procedure described above. The detailed behaviour of the disc near a resonance 
depends, in principle, on a number of factors: the relative importance of various collective effects (such as pressure, self- 
gravity or viscosity), the level of nonlinearity, the vertical structure of a three-dimensional disc, etc. It is known that the 
torque exerted at a Lindblad resonance is independent of these details pro vided that the response is linear and localized near 
the resonant orbit jMever-Vernet fc Sicardv 1987 : Lubow fc Ogilvie 19981 ). while nonlinearity in fact makes little difference 
[ Yuan fc Cassenlll994l ). On the other hand, the torque e xerted at a corotation resonance does depend in an important way 
on the relative importance of nonlinearity and viscosity ( Ogilvie fc Lubow 20031 ). 

We therefore adopt a minimal description of the collective effects in the disc, by introducing viscous behaviour without 
attempting the formidable task of expressing the full Navier-Stokes equation in terms of orbital elements. The dynamical 
equations adopted so far refer to individual particles, and the orbital elements of a given particle are functions of time only. 
To deal with a continuous disc we adopt a semi-Lagrangian description with spatial coordinates (ao,Ao) labelling the fluid 
elements at any instant0 Note that the unperturbed variables satisfy ao = constant and Ao — not = eo = constant for any 
fluid element, with no = (GM/ag) 1 ^ 2 . Therefore (oo,eo) are true Lagrangian coordinates but are less suitable because of the 
rapid shearing in the unperturbed flow. The total time-derivative d/dt translates into the Lagrangian derivative 

D - ■ - (65) 



Dt 



d d 
dt oXo 



Even a simplified model of collective effects has certain minimal requirements. The terms added to the dynamical equations 
should regularize the resonant singularity that occurs in the first-order linearized response. At the same time they should have 
the correct form to conserve angular momentum while dissipating energy. Since viscosity is known to be required to resolve 
the singularity in the linearized equations at a corotation resonance, we add diffusive terms to each equation for the purposes 
of regularization. The diffusion coefficient v can be identified as the effective viscosity of the disc. Specifically, we adopt the 
model 



DA _ dH /dmy 1 d I dm dX \ 
Dt dA \ dao J dag \ dao dao J 



DA 



dH 
' dX 



+ (?) 
Vdao / 



d 
da 



t dm dA \ 
\ dao dao J 

Dz a .dH /dm\ _1 d I dm dz a \ 
Dt dza \ dao J dao \ dao dao J 



) - 2v 


dzp 




da 



(66) 



(67) 



(68) 



7 For a three-dimensional description a third coordinate is required to describe the distance from the 'midplane' of the disc. 
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where m(ao) is the cumulative mass function of the disc, i.e. the total mass interior to the orbit labelled by ao. (For a circular 
disc, dm/dao = 27rEoo, where E is the surface density.) The reason for writing the diffusive terms in this way, rather than 
simply as v(d 2 A/dcio), etc., is so that they properly conserve the mass- weighted integrated quantities. Now if a certain quantity 
obeys a diffusion equation, the square of that quantity is both diffused and dissipated@ By allowing z a to diffuse we therefore 
allow the specific AMD |z/3| 2 to diffuse and dissipate. The same dissipative term is added to the 'energy' equation (|67[) so that 
the specific angular momentum L z = A — \zp\ 2 diffuses conservatively and does not dissipate. 

To see how this model operates, consider the behaviour in the absence of a perturber. With suitable boundary conditions 
we have, after an integration by parts, 



d_ 

dt 



dm = -2 



dao 



dm, 



(69) 



showing the dissipation of the total AMD@ Energy (or A) is also dissipated, but the total angular momentum perpendicular 
to the reference plane is conserved: 



dm = 0. 



In these integrals the mass element is 

dm = -7-4— dao dAo- 

2-7T dao 



(70) 



(71) 



We may neglect the slow effects of diffusion on the smooth unperturbed state. Consider now a first-order equation such 



DAi 



= -Re\ikRe i *°] + (— ) 

L - 1 V da / da \ da da J 



(72) 



Unlike the ordinary differential equation (|12|l . this is now a partial differential equation in variables (ao,Ao,i). Since the 
equation is linear and the forcing term depends on Ao and t only through the exponential, the solution is of the form 
Ai = Re [Ai(ao) e"* ] , where Ai satisfies the ordinary differential equation 



i^oAi = ikR + 



/dm\ _1 d / dm dAi 



V dao / dao V dao dao 



(73) 



It is clear that the viscous term is required to regularize the solution at the resonance where (j>o — 0. For a localized response 
the variation of all coefficients with ao can be neglected, except for (j>o ~ Dx, where x — ao — ag is the distance from the 
resonance and D = d(j>o/dao is the rate of detuning. Then Ai = fkR, where f{x) now denotes the solution (decaying as 
x — > ±oo) of the rescaled equation with unit forcing, 



d 2 / 

Standard Fourier-transform methods (cf. iMever-Vernet fc Sicardvlll987r i give 



(74) 



fix) 



i 



H(-ksgnD) dk, 



(75) 



where H is the Heaviside step function. Note that / ; 
for all x (Fig. 1) if s » (uD 2 ) 1/3 . Furthermore 

2 



(Dx) 



h 1 far from the resonance. In fact / resembles (Dx — is) 



fidx 



dx 



dx 



7T 

W\ 



(76) 



Therefore the resolution of the resonant singularity by viscosity gives a solution with properties very similar to that of the 
simple replacement (|59[) . In particular, the integrated effect of the imaginary part of /, which gives rise to the irreversible 
behaviour, is independent of v and identical to that obtained using equation (|59[1 . 

The solution for Ai and z a i follows in the same way, with the complex f(x) replacing the (fig 1 of the original analysis of 
Section [2]20 The solution proceeds similarly to the second order, with one important exception. The dissipative term in the 



8 Specifically, if dtu = d xx u and v = \u\ 2 , then dtv = d xx v — 2\d x u\ 2 . 

9 A slightly more sophisticated model would allow the dissipation to be related to gradients of z a multiplied by some function of ao, 
e.g. if the gradient of T, rather than that of Z2 , is the important quantity in a warped disc. However, such refinements are inconsequential 
in this context. 

10 It is easy to generalize this treatment of collective effects to allow different viscosities to act on the mean-motion variables (A, A), the 
eccentricity variable zi and the inclination variable 22- Complex viscosities are also permissible, giving solutions in the form of attenuated 
waves emitted from the resonance. 
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Figure 1. Real and imaginary parts of the complex function f(x), which represents how the resonant singularity at x = is resolved. 
The solid line shows the viscous solution, equation (175 D . while the dashed line gives the simpler 'damped' model / = (Dx — is) - 1 for 
s = 1. Without loss of generality, units are chosen such that D = v = 1. 



'energy' equation (|67[) provides an additional contribution to (A), which is needed to enforce angular momentum conservation. 
Using equation (|76|l and carrying out an integration by parts, we find 



(A)d dm = — — 

x ' 2D 



O + i) 



OR 



— -(k 
dao 



dR 



dm ,2, D |2 d / dm 



da 



+ k z \R\ 



dao VdAo 



(77) 



Here (A)d includes the dissipative contribution, which has the effect of changing k into k ± 1. If the disc remains circular, 
A is the specific angular momentum and the above rate of change of actio n corresponds to the torque on the disc. It agrees 
in form with the torque formulae for Lindblad and corotation resonances ijGoldreich fc Tremaine! fl979). Note that dm/dAo 
( = 47rS/n) is inversely proportional to the vortensity of a circular Keplerian disc, so the corotation torque is proportional to 
the gradient of inverse vortensity while the Lindblad/vertical torques are proportional to the surface density. 

It appears that th i s pro cedur e cannot be followed if action-angle ecc entricity and inclination variables are used as in 
Goldreich fc Tremaine ( 1981 ) and Borderies. Goldreich fc Tremaine] 1 1984 ). We have found it essential to allow eccentricity 
and inclination to diffuse (or propagate) at first order, to resolve the resonant singularity, and for the accompanying dissipation 
to appear at second order. This is most naturally achieved by using complex variables for eccentricity and inclination, to which 
the energy and AMD are related quadratically. 



3 THE UNRESTRICTED PROBLEM 
3.1 Canonical equations 

We now extend the problem to treat dynamically the orbit of the companion of mass M' around the central mass M. We 
adopt canonical variables (A', A', q' a ,p' a ) for the two-body problem, where A' = [G(M + M')a'] 1//2 and q' a and p' a are defined 
by analogy with equations (|35[l - (|38p . 

The Hamiltonian of the unperturbed two-body problem is ^Hq, where \i = MM' /(M + M') is the reduced mass and 



H'o 



[G(M + M')f 



G(M + M') 



2A' 2 2a' 
giving rise to an unperturbed mean motion 



\ ' > r \>\ dH 



G{M + M') 



1/2 



(78) 



(79) 



The total Hamiltonian of the coupled system is then of the form 
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H = J {Ho(A)-eRe[R{A,q ct ,p ol ,A',q' a ,p' a )e l<1 ']} dm + fiH' (A') 



and Hamilton's equations read 
A = n(A) — eRe 
A = eRe [ifeiJe 10 ] , 



OR h, 

P. 

dA 



q a — —eRe 



dR cH 

dp a 



eRe 



dR c i 4 
dq a 



as for the restricted problem, together with 



dR i4 
OA' 6 



A' = n'(A') - - / Re 
A 4 



A' = ~ J Re [ik'Re**] dm, 



dm, 



Re 



p' a = - Re 



dR ci0 

dp' a 

dR i4 , 
dq' a 



dm, 



dm, 



(80) 

(81) 
(82) 

(83) 
(84) 

(85) 
(86) 
(87) 
(88) 



determining the feedback on the two-body orbit. 



3.2 Perturbative solution 

The solution can be developed in powers of e as in Section T2.2I When it comes to calculating the secular evolution at 0(e 2 ) 
there are no cross-contributions from different regions of the disc because their first-order oscillations are out of phase. We 
then find that the secular dynamics at second order, away from resonance, is governed by the Hamiltonian 



Ti. dm 



(89) 



where Ti is given by equation (I31|l or equation (145[l . as in the restricted problem, except that TC is now a function of 
(A,q a ,p a ,A',q' a ,p' a ). In detail, 

dH 



dH 

dZa 



1 1 /\ 1 f dH , 

1 f dH 
i(«a) = - / "q r" dm. 
fj, J dzj 

The total AMD of the system is conserved under this precessional dynamics because 

dH 



( / |zs| 2 dm + n\z' p > 2 



*0 



dzp 



„dH , dH dH . , 



■/9 



3^ 



dm = 0, 



(90) 
(91) 
(92) 
(93) 

(94) 



where 9 is an angle of rotation of the entire system about the a-axis. (Under a rotation through 59, vj and Q decrease by 89 
and so 8z a — —\z a 89.) This conservation law is in addition to the conservation of J Hdm itself. 

Resolution of the resonant singularity can again be achieved using the replacement (|59[l or a more sophisticated version 
thereof. We find 



(A) = -ifc/i 





dR 


2 


dR 






dzp 




dz* 





(95) 
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d 2 R* OR 



d 2 R* dR 



dzadz*p dzp dzadzp dz*p J 4 ydzadzp dzp dzadz*p dzp 



d 2 R dR* 



d 2 R OR" 



4 k dX [ fR d* +f R d*)' (96) 



2 ■" 



dR 



dzp 



dR 



dm, 



(97) 



'OA 



a 2 ir ag _ 

dz'adzp dzp 



d 2 R* dR 
dzadzp dzi 

, dR 
dz'* 



dm. 



d 2 R dR* 



d 2 R dR* 



dzadzp dzp dzadzp dzp 



(98) 



Equation (|95|l requires modification to account for dissipation, as discussed in Section \2. 61 leading again to equation (|77|) . 

According to this analysis the secular behaviour of the disc is governed by the precisely same equations as in the restricted 
problem, but the orbital elements of the companion now also evolve in a consistent way through the feedback of the disc 
on the two-body orbit. This coupled dynamics is of special interest when the disc and companion have comparable angular 
momenta, as may occur in protoplanetary systems or planetary rings. 



4 APPLICATIONS 

We now work out some of the implications of the preceding analysis. We restrict attention initially to effects of first degree in 
eccentricity and inclination, and then consider an example of nonlinear behaviour in Section \4. 41 



4.1 Second-order resonances 



In Section 12.51 we considered the behaviour of eccentricity at an interior second-order resonance in the circular restricted 
problem. We now generalize the result to permit both the disc and the companion to have small eccentricities and inclinations, 
and to allow dynamical feedback on the companion's orbit. 

For an interior second-order j : j — 2 resonance we have, to lowest degree in eccentricity and inclination, 



i<ii) = 7 



1 / GM' \ 2 / 2 \ 3 / 



V a ' 
1 f GM' 



VA 



— /[2/ 45 (2/ 45 f + /«£')] . 



4 ^ 



K4) 



GM' 



GM' 



1/2 



1/2 



Ac 



/[/m(2/«£ + /ibO] dm, 



dm. 



(99) 
(100) 
(101) 
(102) 



These expressions, which follow from the terms of second degree in the disturbing function (| A5|) . derive only from the 
'Lindblad/vertical' aspects of the resonance; there are no 'corotation' terms here. The rate of change of total AMD is 

1 / GM' x 



\zp\ 2 dm + n\z'p\ 2 



2/i( \2f 45 £ + f 49 £'\ 2 + -f 2 7 \I-l'\ 2 



dm. 



4\ a> 

For a localized resonance /i is strongly peaked and we may use equation (|76[1 to express the result as 



GM' 



n dm 
AolDI da 



\2f 45 £ + f 49 £'\ 2 + ^ 7 \l-l'\ 2 



(103) 



(104) 



evaluated at the resonance. The equivalent result for an exterior second-order resonance is 



7r dm 
AolDI do 



|2(/s 



\a- 2 5j, 3 )£ + f 49 £'\ 2 + -f 2 7 \I-l' 



(105) 



/ GM' 

V a 

The meaning of these results is that second-order resonances allow an exp onential growt h of AMD for small eccentricities 
and inclinations. We have already explained the connection with the result of iLubowl (|199lh concerning the growth of eccen- 
tricity in the disc when the companion has a fixed circular orbit. (It is noteworthy that the indirect term affects, and largely 
cancels, the eccentricity growth rate associated with the 1 : 3 resonance in an exterior disc, so it should always be included 
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in studies of satellite-disc interactions.) On the other hand, if the disc is fixed artificially to be circular, the eccentricity of 
the companion can grow exponentially. Using the fact that f4g « — (j 2 /27r)[5A"o(4/3) + (19/4) Ai (4/3)] for j ^> 1, where K is 
the modified Bessel f unction, we find a growth rate for £ ' that agrees with equation (8) of the analysis of eccentric Lindblad 
resonances bv lWardl l|l988N "l in the case of a circular disc and e'<l. 

For the first time, though, we see here the effect of the resonance in the general case when both the disc and companion 
are eccentric. Certain linear combinations of the complex variables £ and £ ' are involved, which depend on j. The case of 
inclination is easier to understand, as the releva nt combination is always X — X' , the complex 'mutual inclination' of the disc 
and compa nion. Indeed, Lubow fc Qgilviel ( 2001 ) we re able to write down such equations for inclination dynamics based on the 
formulae of Borderies. Goldreich fc TremaineT i 1984 ) (for a fixed disc) and the idea that only the c omplex mutual inc li natio n 
could enter. Using the fact that f$7 — (a/2)b^Z 1 \a), we find agreement with equation (48) of Lubow fc Ogilvie ( 2001 ). 
Although 2/4 5 //49 7^ —1, it approaches —1 for large j, and so a similar combination £ — £' is involved in that limit. Especially 
in a protoplanetary system in which the disc and planet (s) have comparable orbital angular momenta and engage in coupled 
secular dynamics, the eccentricities and inclinations of both disc and companion(s) must be considered simultaneously. 



4.2 First-order resonances 



In the case of a circular and coplanar system the only operative term in the disturbing function for an interior first-order 
j : j — 1 resonance is 



R=^fhr£. 

n ' 



(106) 



This leads to no evolution of z a or z' a (which remain zero) but does produce a torque on the disc, which is (according to 
equation [77)l 

2 

(107) 



Aol-Dl 



GM' 



f 2 dm 
da 



Using the fact that f 2 7 = — jb^ 2 (a) — (a/2)db^ 2 /da, this torque can be shown to agree with equation (13) of Goldreich fc Tremainel 
1 19801 ). The equal and opposite torque on the companion follows from equation (|97[l . 

When small eccentricities and inclinations are allowed for we obtain the following dynamics to lowest degree: 
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1/2 
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Mo 



7 f(f 31 -2a5 jl2 ) 



l dm, 
(110) 

dm. (Ill) 



Both 'Lindblad' and 'corotation' aspects of the resonance enter here. Regarding the 'corotation' terms (those involving 
derivatives with respect to Ao), we again have precessional behaviour that diverges, this time more strongly, at the resonance, 
owing to the appearance of d//dAo. Resolution of the resonant singularity leads to an eccentricity growth rate (e.g. when 
£' — 0) that changes sign across the location of the resonance. Which sign is obtained for the disc as a whole depends on 
whether dm/dAo is greater on one side or the other. The net effect for a narrow resonance is proportional to d 2 m/dAo, or 
equivalently d(E/n)/dao, as is characteristic of corotation resonances. 



1 In his equation the second K\ should be Kq. 
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It is notable that (22} and (i 2 ) do not vanish when the disc and companion have inclined but coplanar orbits (X = X' 7^ 0). 
Using the relation fa — fa = 2^/27 it is possible to show that (22) and (z' 2 ) in this case can be attributed entirely to (A) and 
(A') while I and I' remain constant. In general, the rate of change of total AMD is 

GM' V (j - 1)tt d /dm\ 



2\D\ dao VdAo/ 



fa£ + (/31 - 2a5 ja )£ \ . (112) 



Note that only the 'corotation' terms contribute to this expression. The equivalent result for an exterior first-order resonance 
is 



GM' 



2 



2\D\ da VdAo/ 



W + (/si - \a-*5^)£\\ (113) 



Using the fact that /31 « —fa ~ (j/7r)[2Jfo(2/3) + Jf 1 (2/3)] for j 3> 1, we f i nd a growth (or decay) rate for £' that agrees 
with equation (12) of the analysis of eccentric corotation resonances by I Ward! (|1988T ) in the case of a circular disc and e' <C 1. 
Whether the AMD grows or decays depends on the sign of the vortensity gradient. In general, as for second-order resonances, 
the eccentricities of both disc and companion must be considered simultaneously, as the resonant effect depends on a linear 
combination of £ and £' , in fact £ — £' in the limit of large j. 



4.3 Zeroth-order resonances 



Zeroth-order resonances are coorbital and are active only if a clean gap is not opened by the companion's orbit. In the case 
of a circular and coplanar system the only operative term in the disturbing function for a zeroth-order resonance is 



CM' 



(114) 



This leads to no evolution of z a or z' a (which remain zero) but does produce a torque on the disc, which is (according to 
equation I77[l 

'gm ,x 2 



itj 2 GM'\ 2 d fdm\ 

m\(~ J (2/i ^ i} d^ldAj 



(115) 



Using the fact that /1 = %b^) a {a), this torque can be shown to agree with equation (14) of Goldreich fc Tremaine ( 1980T ). and 
the equal and opposite torque on the companion follows from equation (|97|) . Setting a = 1 renders the coefficients singular, 
and this must be resolved in practice by a smoothing process that represents an averaging of the potential over the vertical 
extent of the disc, which we do not consider here. 

When small eccentricities and inclinations are allowed for we obtain the following dynamics to lowest degree: 



i<ii) 



' " 1 (ij (2/2 + ^ a5j ' l] [ (2/2 + 2 a ^M)£ + iho - aS j , a )£ r ] + f [(2/2 + \a5 jA )£ + f 10 £'] } , (116) 



V a ' 



1 / GM' 
4 ^ a' 



i<4> = 



GM' 



1/2 



Ai j (ao") { " f[ ~ ho ~ aSj > 2) [ (2/2 + 5 a5 ^ £ + - a5 ^) £ '] 
+ffa [(2/2 + \a5 jtl )£ + faS 1 ] } dm, 

1 1 i (^") 2 (^) 1/2 © 1 " mh 2aSj - l] [i(2/3 + a5j ' l)J+ K/l4 " 2a5ja):r ' ] 

+/' |/i4 [|(2/3 + c^,i)X+ i/i 4 l'] } dm. 
For brevity we omit 'corotation' effects here and focus on the 'Lindblad' terms. The rate of change of total AMD is 
_ / GM' 



(117) 



(118) 



(119) 



7T dm 



|(2/ 2 + \a&i,i)e + (fa - aSj, 2 )£'\ 2 + |(2/ 2 + \a5 Jtl )£ + fa£' 



Ao\D\ da Q 

+\\{2f 3 + a5 j , 1 )T+\{fa-2a8 j , 1 )T'\ 2 + |±(2/ 3 + aS 3<1 )l + \fal'\ 



(120) 



again omitting 'corotation' terms. Althou g h the se expressions cannot be consistently evaluated without a softening procedure 
and the further considerations of Ward ( 19881 ). we see that the effect of the coorbital 'Lindblad' terms is to damp the 
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eccentricity and inclination of the disc or companion if the other is fixed at zero eccentricity and inclination. Generally, 
however, the effect of the resonance depends on linear combinations of £ and £ ', and of I and T' , which are again £ — £' 
and I — T' in the limit of large j. In this sense it is the (complex) differences in eccentricities and inclinations of the disc and 
companion that are being damped, so these terms attempt to equalize both e and vj (or i and SI) of the disc and companion. 

Again, it is notable that (22) and {z' 2 ) do not vanish when the disc and companion have inclined but coplanar orbits 
(X = T 7^ 0). Using the relations 2/3 + /14 + a8j,i = —2/3 — /14 + ct5j,i = —j(2fi — a<5j,i) it is possible to show that (22) and 
{z'2} in this case can be attributed entirely to (A) and (A'} while T and T' remain constant. 



4.4 Effects beyond the first degree in eccentricity 

For a companion on a circular orbit (e.g. in a close binary star) and a coplanar disc, the disturbing function for an interior 
second-order resonance is 
GM' 



R 



^[U,£ 2 + {Us~\U.)\£\ 2 £ 2 + 0{£ 6 )] 



[121) 



In this case 



4 y 

+0(z°). 

The effect on the AMD of the disc is 



I |2 
Zl Zl 



(122) 



\zfj\ dm 



+0(z b ), 



( GM' 



dm 



MDl 1 d> t 4 ^! 2 + U <*« - ^l £ | 4 ] + « - 2)Ao^ (£) /&|fi|< 



(123) 



evaluated at the resonance. 

As a specific example, for the 3:1 resonance we have /4s = 0.5988 and fee 
equation are therefore 

. dm 



1.434e 2 ^ [1 - 1.470 C 2 - \pe 2 + 0(e 4 )l , 
dao 1 - 



-0.1936. The terms in braces in the above 



(124) 



where /3 is the logarithmic derivative of vortensity with respect to a. This result implies that the eccentricity growth rate 
is substantially reduced when e is not small, providing a possible means to limit the growth of eccentricity in superhump 
systems. 

To evaluate consistently the evolution to higher degrees in eccentricity and inclination, further resonances must be taken 
into account. For an nth-order resonance (n > 2) we have R oc £ n , leading to a growth of eccentricity in the disc with 

£ oc \£\ 2n - 4 £. 



5 SUMMARY AND CONCLUSIONS 

We have presented a new and general analysis of mean-motion resonances between a Keplerian disc and an orbiting companion. 
The emphasis of this treatment is to provide a systematic method to calculate the rates of change of eccentricity and inclination 
variables of the disc and companion associated with resonances of various orders. 

Near a mean-motion resonance, terms in the disturbing function corresponding to the appropriate commensurability 
produce an oscillation of the orbital elements of the disc and companion at first order in the perturbation. At second order, 
a secular effect is obtained in the form of a precessional dynamics that diverges as the resonance is approached. When the 
divergence is resolved by taking into account the (dissipative) collective effects of the disc, a finite precession rate is obtained, 
together with an irreversible evolution, such as a growth or decay of eccentricity or inclination. 

Our work differs from earlier methods in which the perturbing potential is decomposed into rigidly r o tating comp onents 



Uur work diners trom earlier methods in which the perturbing potential is decomposed into rigidly rotating components 
and the eccentricity evolution is deduced indirectly from resonant torques ( Goldreich fc Tremaine 1980l ; Ward 1988h . Such 



methods do not apply to general satellite-disc interactions in which both the disc and companion have eccentric and/or 
inclined orbits. Where appropriate, however, we obtain full agreement with this earlier work. By using the classical disturbing 
function, we benefit from the ready availability of expansions to high degrees in eccentricity and inclination. 

After considerable thought, we recommend the use of complex canonical Poincare variables, as defined in Section 12.31 
to describe the eccentricity and inclination dynamics. These allow the most compact expression of the dynamical equations 
and the disturbing function. They are also directly related to the angular momentum deficit (AMD), which is a quantity of 
fundamental importance, being a positive-definite measure of the departure of the system from circularity and coplanarity 
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that is conserved in first-order secular interactions. By concentrating on the evolution of the system's AMD we can separate 
the effects that lead to irreversible growth or decay of eccentricity and inclination from those that are merely precessional 
in character. Moreover, the use of complex variables leads, in the case of small eccentricities and inclinations, to linear 
evolutionary equations that would be extremely cumbersome in any other representation. 

This work does not address certain refinements that are necessary for a more complete description of resonant interactions. 
The important omitted effects are nonlinearity (in particular the saturation of corotation resonances), torque cutoff effects, 
resonant shifts and vertical averaging. It is likely that all of these could be addressed within the present framework. 

The equations derived here do not stand alone but should be combined, in future work, with the complex one-dimensional 
partial differential equations describing the s ecular evolution of eccentricity and inclination in t he disc. In simple cases such 
a programme has already been carried out ( Lubow fc Ogilviel 2001 ; Goodchild fc Ogilvie 2006h . For small eccentricity and 
inclination the outcome is a set of secular normal modes, being a continuum analogue of the Laplace-Lagrange secular theory, 
which describe rigidly precessing configurations of disc(s) and companion(s). The contribution of mean-motion resonances 
allows these modes to grow (or decay), and any growth m ust be offse t again st the viscous damping of the mode within the 
disc. Coupled eccentric modes have also been discussed by iPapaloizou l|2002h . 

Previous discus sions of eccentricity gr o wth i n protoplaneta ry s ystems, for examp l e, hav e presented an interesting but 
incomplete picture. Goldreich fc Tremainel ( 19801 ). Ward! 1 19881 ) and Goldreich fc Sari (2003) discuss the close competition 
between eccentric Lindblad resonances and eccentric corotation resonances. This is equivalent to comparing equations (|104[) 
and (|112|) with the complex eccentricity £ of the disc set to zero (and X = X' = 0). Since the disc becomes eccentric through 
secular (and possibly resonant) interactions, the relative weighting of all these resonances is adjusted in a way that depends 
on the eccentricit y distr ibution within the disc. A self-consistent solution of the coupled secular dynamics, as carried out by 
Lubow fc Ogilviel ( 2001 ) in the case of inclination, is therefore required to address the issue. It rem ains to be seen whethe r 
protoplanetary systems support growing eccentric modes in the presence of viscous damping (see also Latter fc Ogilviel 20061 ). 

Another application of this work is to close binary stars, where accretion discs can someti mes become e ccentric through 
the action of the 3:1 resonance. This work suggests various elaborations of the description by iLubowl (|199ll ). In addition to 
eccentricity growth, the 3:1 resonance provides an second-order precession that may affect the shape of the eccentric mode and 
alter the global precession rate of the disc. The growth rate is weakened when the eccentricity is no longer small, leading to a 
possible saturation mechanism for the superhump instability. Finally, higher-order resonances such as 4 : 1 provide eccentricity 
growth of the form £ oc \£\ 2 £ which, although not leading to a linear instability, might conceivably support a global eccentric 
mode in binary stars of larger mass ratio, where the disc is too small to contain the 3:1 resonance. 
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APPENDIX A: EXPANSION OF THE DISTURBING FUNCTION 



In this section the disturbing function is given for various commensurabilities, correc t to fourth degree in eccentricity and 
inclination. The coefficients fa are defined in Appendix B of lMurrav fc Dermott 1 1999h . 



Al Secular terms 



|2T - X'| 2 + | (\£\ 2 + \£'\ 2 ) ill" +1*1') - | (\£\ 2 \1\ 2 + |£'| 2 |T'| 2 )] 



*o.o = ^{/i + / 2 [|£| 2 + |£'| 2 

+(A - i/2)|£| 4 + A [\£\ 2 \£'\ 2 - (1^ I 2 + |£f) \1-1'\ 2 ] + (A - lf2)\£'\ 4 + \fn (£ 2 £" 2 +£* 2 £' 2 ) 
+i (££"+£*£') [ho + (fn - |/io)|£| 2 + (/i2 - §/io)|£f + \hz\l-lf] + ±f a \l-?\ 4 
+^(/ 8 -4/ 5 ) [(\1\ 2 + \1'\ 2 ) {11" +T1') -A\l\ 2 \l'\ 2 } + |/ 18 [£ 2 (T — l") 2 + £* 2 (1 — 1') 2 ~\ 
4-g/io [££'(T - I") 2 +£*£"(! -I') 2 ] + i/ 20 [£' 2 (T - 1") 2 + £" 2 (1 ~ I') 2 ] 

-5(2/19 + As) (££"-£*£') (11" -1*1') }, (Al) 

using coefficients f j as defined in Tabl e B.l of lMurrav fc Dermotj dl999l ) with j = 0. This expression can be shown to agree 
with that given bv lLaskar fc Robutel jl995l ). Note that Ro,o is real, and contributes negatively to the secular Hamiltonian. 



A2 Zeroth-order resonances 

Co-orbital j : j resonance with j > 1, for which k = j , k! = — j , <^o = jn — jn' and a = 1: 

R = ^r-{(2/i ~ «*i.O + (2A + |««5 fll ) (|£| 2 + |£'| 2 ) + (2/ 4 - i/ a - £a* w )|£| 4 + (2/s - |a<5,M)|f \ 2 \£'\ 2 
+(2/ a - i/ 2 - ^a^Olf '| 4 + /i 7 £ 2 £:'* 2 + (/it - &a*j,i - §a^, 3 )f * 2 f ' 2 
+££" [Ao + (Ai - gAo)|£| 2 + (/12 - |Ao)|f'| 2 + |/i3 (|X| 2 + |T'| 2 ) + |/ 2 2T*I' + ifeXX'*] 
+fi*fi'[(/io - aSjj,) + (Ai - |Ao + |a<Si,2)|f I 2 + (/12 - |Ao + Iol8 ];2 )\£'\ 2 + + a^, 2 ) (|X| 2 + |X'| 2 ) 

+ i/ 2 2XX'* + K/23 - 2a5 J , 2 )J*J'] + 1(2/3+0^,1) (|X| 2 + |X'| 2 ) + i/ 8 (|X| 4 + |J'| 4 ) + i(2/„ - a<5 3 ,i)|I| 2 |X'| 2 
+i/ 26 T 2 T'* 2 + i(/ 26 - aS^T 2 !' 2 + \11" [/14 + (/i 5 + i/14) (|f I 2 + |£f) + |/i6 (|X| 2 + |X'| 2 )] 

J' [(/u - 2a^,r) + (As + |/i4 + §a* ifl ) (|£| 2 + |£'| 2 ) + i(/ M +a6 jtl ) (\1\ 2 + \l'\ 2 )] 
+±(A + 2/ 7 ) (j£| 2 |T| 2 + |ff |Tf ) + I(2/ 7 - 1^,0 (|£| 2 |J'| 2 + |£'| 2 |X| 2 ) 
+±(X 2 +X' 2 ) (/igf* 2 + fi 9 £*£'* + f2o£" 2 ) + \{T 2 +1" 2 ) (fis£ 2 + fytSS' + ho£' 2 ) 
+\11' (f2i£* 2 + f24£*£" + f2 5 £" 2 ) + \T1" (f~2i£ 2 + f24££' + f 25 £' 2 ) 

-±0163,1 [£* 2 {1- 1') 2 +£' 2 (T -T'*) 2 ] }. (A2) 

Here / denotes the coefficient /, in which j is replaced with — j. Setting a = 1 renders the coefficients singular, and this must 
be resolved in practice by a smoothing process that represents averaging over the vertical extent of the disc. 

A3 First-order resonances 

Interior j : j — 1 : resonance with j > 2, for which k = j — 1, k' = —j, 4>o = (j — l)n — jn' and a — a/a' = [(j — l)/i] 2,/3 : 

R = ^f{f2r£ + (Ai - 2a5j, 2 )£' + (/ 28 - \f„)\E \ 2 £ + /»|£'| B £ + (A2 + aS j<2 )\£ \ 2 £ ' + (/ 33 - |Ai + |a^, 2 )|£f £' 
+f 35 £ 2 £'* + (Ae - f oc&^ z )£' 2 £ * + i[f 39 £ + (/ 42 - 4a^, 2 )f] [T {i -1) -l'(l" -T)} 

+Uho£ + /43O [I(I'*-I*)-I"*(T'-T)] +i(/3 7 f* + A8f'*)(2:-2:') 2 }- (A3) 
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Exterior j — 1 : j resonance with j > 2, for which k = — j, k' = j — 1, <j>o = —jn + (j — l)n' and a = a' /a = [(j — l)/i] 2 ^ 3 : 
CM' 

R = —{hr£ ' + (/si - hor 2 S jfl )S + (/m - |/ 27 )|£'| 2 £' + / 29 |£ | a f + (/s2 + |a" 2 ^, 2 )|£ '| 2 £ 
+(/33 - I/ai + ^a- 2 S jt2 )\£\ 2 £ + f 3o £' 2 £* + (/ 36 - \a- 2 5 j>3 )£ 2 £'* 
+\\hs£' + (h2-a- 2 6j,2)£] [l'*(l-l')-l(l* -I'*)] 

+g(/4o£' + /«£) [2* (I* —l'*) — X* (Z — Z')] + i(/s7f" + /ssfi *)(Z' -X) 2 }- (A4) 
A4 Second-order resonances 

Interior j : j — 2 resonance with j > 3, for which = j — 2, fc' = — j, 0o = (j — 2)n — jn' and a = a/a' = [(j — 2)/j] 2 ^ 3 : 
CM' 

R = —{Ur,£ 2 + f4 9 ££' + (fa - f a5j, 3 )£' 2 + (he - ±/«)|£ | 2 £ 2 + /4t|£ '| 2 f 2 + (/so - | 2 £f 

+(/si - |/4 9 )|£'| 2 ££' + (/ M + ^aSj, 3 )\e \ 2 £ 12 + (/ M - i/53 + fa^iefg 12 + / 68 £ 3 £'* + (/ 69 - *a6 jA )£ *£ 13 

+\hj(i-i'f + ±f 8 o[i 3 (i" -T) +i'*(i 3 —i' 3 )] + i/ 8 i[T' 3 (r — 1'*) +i*(i ,3 -i 3 )] 

+-L/ 6 i[X 2 X'(X'* -X*) - 3XjX'| 2 (X' -X)] + i/ 67 [X' 2 X(X* -X'*) - 3X'|X| 2 (X - X')] 
+ ±/ 57 (X-X)(|£| 2 X- |£'| 2 X) + I[/ 72 £ 2 + f 74 ££' + (frs - f a5 3 , 3 )f' 2 ] [X* (X - X) - X (X* -X*)] 
+g(As£ 2 + /75^' + /79f ' 2 ) [X(X* -X*) -X'*(X -X)] + i(/ 58 |£| 2 + / 59 |£'| 2 + Aoff + fri££'l(l-l') 2 }. (A5) 
Exterior j — 2 : j resonance with j > 3, for which k = — j, k! = j — 2, <j>o = —jn + (j — 2)n' and a = a' /a = [(j — 2)/j] 2 ^ 3 : 

R = ^{Ub£' 2 + f m e'£ + (h, - %a- 2 8j, 3 )£ 2 + (he - \hs)\£'\ 2 £' 2 + Ui\£\ 2 £ a + (/so - IM\£'\ 2 £'£ 

+(/ M + ±cr 2 8 jfl )\£'\*£* + (ho - i/53 + ^a- 2 8 3 , z )\£\ 2 £ 2 + / 68 £ ,3 £* + (/ 69 - l^ 2 8 ]A )£'* £ 3 

+±/ 57 (X -X) 2 + i/ 80 [X ,3 (X* -1") + T(l' 3 ~l 3 )] + i/ 81 [X 3 (X'* -T) +1'*(1 3 -I' 3 )] 

+±fei[l' 2 l(T -1") - 3X'|X| 2 (X-X')] + -^/ 67 [X 2 X'(X'* -X*) - 3X|X'j 2 (X -X)] 

+±hv(l' -1)(\£'\ 2 1' - \£\ 2 1) + WhzE' 2 + f 7i £'£ + (Jn - ^S.^S 2 } [l"> (1 -l') -1(T -I")] 

+Hf™£' 2 + frs£'£ + fr S £ 2 ) [l'(T -!'*)-!* (1-1')] + i(hs\£'\ 2 + h & \£\ 2 + fn£"£ + fnS'S *)(!' -I) 2 }. (A6) 

A5 Third-order resonances 

Interior j : j — 3 resonance with j > 4, for which k — j — 3, k! = — j, <j>o = (j — 3)n — jn' and a = a/a' = [(j — 3)/j] 2 ^ 3 : 
CM' 

R = ^7- [ftae 3 + /asfi 2 f ' + /84£ £ ' 2 + (/ 85 - faS jA )£' 3 + \(f m £ + hr£')(l - X') 2 ] . (A7) 

Exterior j — 3 : j resonance with j > 4, for which k = —j, k! = j — 3, cj>o = — jn + (j — 3)n' and a = a' /a = [(j — 3)/j] 2 ^ 3 : 
CM' 

R = ^- [fta£ ' 3 + hi£' 2 £ + h4£'£ 2 + (hs - \a~ 2 5 jA )£ 3 + l(f 86 £' + f S7 £ )(X' - X) 2 ] . (A8) 



